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, The impact of nonrenormalizable gluon operators upon inclusive jet cross sections 

is studied. Such operators could arise in an effective strong interaction Lagrangian from 
gluon substructure and would induce observable cross section deviations from pure QCD 
at high transverse jet energies. Comparison of the theoretical predictions with recent CDF 
data yields a lower limit on the gluon compositeness scale A. We find A > 2.03 TeV at 
95% CL. 
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The inclusive jet cross section data from the 1988-89 Fermilab Tevatron run span 
seven orders of magnitude and include the highest transverse jet energy measurements 
reported to date [Q. These data provide a stringent test of quantum chromodynamics 
and constrain possible new physics beyond the Standard Model. In particular, they set 
improved limits on hypothetical quark substructure. At energies small compared to the 
compositeness scale A, the dominant effects from quark substructure can be reproduced by 
four-quark operators in a low energy effective Lagrangian 0. Although their coefficients 
are unknown in the absence of a detailed theory of preon dynamics, the general impact 
of these nonrenormalizable operators upon parton scattering may be estimated. The re- 
markable agreement between the experimental measurements and the predictions of QCD 
then places a bound on the quark compositeness scale. CDF finds a lower limit on A of 
1.4 TeV at the 95% confidence level ||. 

In this letter, we reinterpret the CDF data to probe for signals of new physics that 
could arise in the gluon sector. Specifically, we consider the impact upon the inclusive 
cross section measurements of nonrenormalizable gluonic operators which may appear in 
the effective strong interaction Lagrangian. Such operators could originate from a number 
of different sources. For example, suppose there exist new heavy colored bosons or fermions 
beyond those in the Standard Model. Such particles would induce nonlocal interactions 
among gluons through loop diagrams. The leading behavior of these graphs can readily 
be extracted and reexpressed via an operator product expansion in terms of local but non- 
renormalizable gluon operators. Alternatively, we might speculate that gluons are bound 
states of some more fundamental preon constituents. Then as in the case of composite 
quarks, preon exchange could generate nonrenormalizable gluon interactions. In the fol- 
lowing, we will adopt a model independent approach and not specify the underlying physics 
whose low energy effects are encoded in the effective Lagrangian. Instead, we simply seek 
to place a limit on its characteristic scale A. 

We first enumerate the lowest dimension gluon operators whose scattering effects 
would be easiest to observe. There exist only two independent operators of mass dimension 
d+2md = 4 — e spacetime dimensions which preserve gauge invariance along with C, P 



and T 0: 
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where 

(2) 

= dPG: - d^GP^ + ii^'^gfa,cG',G^, 

and denotes the renormahzation scale. AU other dimension-((i+2) gluon operators either 
vanish or reduce to combinations of the two in eqn. 

The contributions of Oi and O2 to parton scattering cross sections have been studied 
in refs. [^^. Operator Oi mediates gluon-quark as well as gluon-gluon scattering and 
interferes with the pure QCD amplitudes for these same processes. One might expect 
its impact to become pronounced at higher energies due to the sizable gluon content of 
colliding hadrons at small parton momentum fractions. However, the helicity structure 
of the 0(1/A^) amplitude from Oi for gg gg scattering is orthogonal to that of pure 
QCD. So the two amplitudes do not interfere. Similarly, the 0(1/A^) interference terms 
in the differential cross sections for gg qq and the other partonic processes related by 
crossing vanish in the limit of zero quark mass. Thus the presence of Oi in the low energy 
Lagrangian affects inclusive jet cross sections starting only at 0(1/A^). This surprising 
and rather disappointing result motivates us to search for signals of gluon compositeness 
in non-gluonic channels. 

Fortunately, operator O2 mediates quark-quark scattering. The classical equation of 
motion 

flavors 

relates the S-matrix elements of O2 to those of a color octet four-quark operator: 

— E (^7.T.,)(^7^T„,). (4) 

flavors 

This identification should be understood as a relation among S-matrix elements and not 
as a true operator identity ^j. Notice that the equation of motion (H) produces a factor 
of g^ multiplying the operator in . No such extra powers of the strong interaction cou- 
pling accompany the analogous color singlet operators that enter into quark compositeness 
analyses. At typical Tevatron energies, the numerical value for g^ is somewhat larger than 

^ The operator obtained from Oi by replacing the antisymmetric structure constant fabc with 
the completely symmetric symbol dabc violates charge conjugation and is identically zero. 
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unity. So we anticipate that the effect upon quark scattering from the four-quark oper- 
ator induced by gluon substructure wiU be shghtly larger than that from the four-quark 
operators generated by quark substructure. 

To complete our operator basis so that it closes under renormalization, we include 
three more four-quark operators along with Oi and O2 in the effective Lagrangian: 

5 

CeS = ^Qcn + '^Ci{ll)Oi{lJ.) (5) 



where 
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Since the underlying short-distance physics responsible for generating these nonrenormal- 
izable operators is not known, there is a fair amount of arbitrariness in the values one 
chooses for their dimensionless coefficients and the compositeness scale A. We will follow 
the convention adopted in previous quark substructure studies and define A to be the scale 
where the magnitude of the gluon operator coefficients Ci or C2 equals 4% [HH. We set 
the coefficients of the remaining operators at the A scale to zero. 

To determine the operators' coefficients at energies probed by the Tevatron, we evolve 
their values down from the compositeness scale via the renormalization group equation 

/^/C',(^) = E(7")..C,(^)- (7) 
j 

The anomalous dimension matrix 
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describes the mixing for Uf active quark flavors among all the dimension-((i + 2) operators 
in the effective Lagrangian. The last four rows of 7 quantify the running of the four- 
quark operators in our basis and were determined from a straightforward operator mixing 
computation. The entries in the flrst row on the other hand were extracted from the highly 
nontrivial anomalous dimension calculations reported in ref. by Narison and Tarrach 
and in ref. |0] by Morozov. Notice that the triple gluon field strength operator runs only 
into itself and does not mix with any of the other four-quark operators, i Since Oi has 
virtually no impact upon parton scattering and does not mix into operators whose effects 
can be detected, it is essentially invisible at 0(1/A^). We will consequently be unable to 
place any limit upon its associated scale. 

It is convenient to decompose the anomalous dimension matrix as 

(8) 

= l-^SDS-'+Oig^) 

where the eigenvalues of 7 are contained in the diagonal matrix D while the correspond- 
ing eigenvectors are arranged into the columns of matrix 5". The general solution to the 
coefficients' renormalization group equation 

= J] exp / dg-L^^ C,{A) (9) 



can then be rewritten as 



C^.(^) =E(^"')r/i^g([^] ^')(Sn,,C,{A) (10) 

where b = —11/2 + n//3 is the coefficient in the one-loop QCD beta function (3{g) = 
bg^/Siv'^ [l^]. Integrating /3, we obtain the strong interaction fine structure constant 



1 - b log — 

TT Mz 



^ As shown in Morozov's paper, the anomalous dimension of the CP-even operator is the 
same as that of the CP-odd & operator which gained much notoriety a few years ago S. 



We choose the constant of integration in the integrated beta function to be as{Mz) = 
0.118 ± 0.007 rather than the QCD scale since parton energies at the Tevatron are 



typicaUy large compared to the Z scale. 

We now turn to computing the inclusive jet cross section in pp collisions. Neglecting 
higher order multi-jet events, we start with the two-jet differential cross section 



a 



da 



{AB 2jets) = ^ 2xaXbPT fa/A{xa)fb/B{xb)+{A ^ B if a b) -^{ab cd) 

abed 

(12) 

expressed in terms of the jets' pseudorapidities (771, 772), their common transverse momen- 
tum (pt), and the momentum fractions {xa.Xb) of partons a and b inside hadrons A and 
B. The partons' distribution functions fa/A{xa) and fb/B{xb) are folded together with 
the differential cross section da/di for the elementary scattering process ab cd. i The 
product is then summed over all possible initial and final parton configurations. To con- 
vert the two-jet expression (|T2|) into an inclusive single-jet cross section, we integrate over 
the pseudorapidity range of one jet, average the other over the pseudorapidity interval 
0.1 < {ill < 0.7 visible to the CDF detector, and multiply by two to count the contribu- 
tions of both jets to the inclusive cross section: 

1 d^a _ 2 d^a 

Ar] J dr]dpT Ar] J V-00 dr]idr]2dpT' 

The result may then be compared with the measurements reported by CDF [|l[|. 
The lowest order QCD predictions for the parton cross sections 



-^{ab cd) = — 
dt s 



{ab cd) = -^^{ab cd) (14) 



have been frequently documented in the literature [|l2| , pJ3(] . They conventionally include 
initial state color averaging factors and are written in terms of the partonic invariants s, t 
and u. The QCD formulae for E(a6 — > cd) are modified by the nonrenormalizable operators 
in our effective Lagrangian which induce the 0(1/A^) interference terms tabulated below: 



^ Note that the Mandelstam invariant i refers to the cohiding partons rather than the incident 
hadrons. 
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Eta'^„') = i^ + 5(^±^l)£!±(^^^ + 0(J,-) (15a) 

^ 9 t2 9 VA4/ ^ ^ 

- ^ ^ ) = + 9 JX-2 + ^lAij (15^) 

^ qq) = - I — + 



9V {2 u2 J 27 iu 

^ 9A2 V i ^ u J ^ 9A2 I i ^ u 

^ gg) = - ^- + 



9V {2 s2 J 27 si 

^ 9A2 I £ ^ s ) 9A2 V £ ^ s 

.„,^,,).^(3-|-|-|).o(i,), (15.) 



Here g' denotes a quark not identical in flavor to quark q. We have dropped 0(1/A^) 
terms in these formulae since we are not keeping track of any dimension- ((i + 4) gluon 
operators whose contributions to E(a6 cd) are of the same order. We have also neglected 
all parton masses. At transverse jet energies of a few hundred GeV, this should be a 
good approximation except for the top quark whose mass we assume is rrit = 140 GeV. 
Initial state top mass effects may be safely ignored as the top content of colliding protons 
and antiprotons is negligible. However for processes involving ti production, we replace 
eqns. ( p.56| ) and ( p.5e| ) with the heavy flavor QCD cross sections given in ref. and 
incorporate 0{1/A'^) interference corrections: 
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_ -.4 (£ - m2)2 + {u- mlf + 2m?s 8 (C2 + Cs)^^ + (C2 - Cg)? ^ 2C2m2s 
Hqq - tt) - - j-^ + - 

+ 0{^) (16a) 



s2 24 (m| - t) (m? 



u 



1 {ml — i){m1 — u) — 2mf{mf + i) 1 (m| — t)(mf — u) — 2mf{mf + u) 
^6 (m? - t)2 ^6 (m? - ti)2 

?>{ml — i) {ml — u) + ml{u — i) ?> {mf — i){ml — u) + m1{i — u) 



8 s(m? - 1) 8 ^(mi^ - u) 

--^(3 + 2m?— 3„,?^)+0(^). (166) 

Note that coefficient Ci of tfie triple gluon field strength operator enters into eqn. ( |16fe| ). 
But since it does not appear in any other cross section formula at 0(1/A2), operator Oi 
has almost no perceptible effect upon the inclusive cross section. We consequently ignore 
it from here on. 

Combining eqns. (|T^) - ([T6|), we calculate the single-jet inclusive cross section as 
a function of transverse jet energy E^- We perform the computation using the leading 
order parton distribution functions of Morfin and Tung (MT set SL) and the CTEQ 
collaboration (CTEQ set L) [|16| as well as the next-to- leading order functions of Morfin 
and Tung (MT sets Bl and S), Harriman, Martin, Roberts and Stirling (HMRS set B) 
and the CTEQ collaboration (CTEQ set MS). All of these parton distribution functions 



along with several others are conveniently contained within the PAKPDF package ||T^ . 

Representative results obtained from the MT set SL structure function evaluated 
at the renormalization scale = are compared with the experimental data in 

fig. |l]. The solid curve in the figure illustrates the predictions of pure QCD with no 
nonrenormalizable operator interactions. Following the example of the CDF analysis fll], 
we have multiplied the theoretical predictions by a normalization factor n to align them 
with the data. A fit for this constant performed over the region 80 GeV < < 160 GeV 
where effects from any compositeness operator terms are negligible yields n= 1.35 ±0.01. 
The resulting agreement between the shapes of the QCD and experimental cross section 
values is striking. There is however a slight suggestion of discrepancy at the highest 
measured transverse energies where compositeness operator effects would be expected to 
first show up. We therefore plot in the same figure the differential cross sections obtained 
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after setting A = 2.0 TeV and C2(A) = —4n in our effective Lagrangian. The resulting 
dot-dashed curve in fig. |1| quahtatively appears to fit the data shghtly better. 

To be more quantitative, we perform a least squares fit for the compositeness scale 
A. First we multiply the differential cross section in each CDF bin by the integrated 
luminosity and bin-width to convert into number of events: 



A' = 4200„b-'x(A£./d,-^)„b. (17) 



We then examine the transformed error bars to determine the statistics obeyed by the 
binned events. The statistical uncertainties in each bin with transverse energy < 
115 GeV are significantly greater than We therefore exclude these non-Gaussian 

points from our least squares analysis. Bins in the intermediate energy range 115 GeV < 
Et < 300 GeV contain large numbers of events and follow Gaussian statistics. At the 
highest transverse energies, the data bins have fewer than 25 events and are described 
by Poisson statistics. We assign Gaussian error bars to these last points following the 



discussion in ref. [|T9[ and then treat them like Gaussian bins containing greater numbers 
of events. 

In addition to the statistical fluctuations associated with each bin, we need to consider 
the systematic errors. Normalizing the theoretical number of events by the factor n removes 
a large, correlated systematic uncertainty. We therefore subtract an averaged percentage 
uncertainty from all the bins' error bars. The corrected systematic errors are then small 
for bins with transverse energies above 80 GeV. 

We adopt the function 

X' = J2^^(V-%A, (18) 

where = A^*"^ — A^^'^'^ represents the difference between the theoretically expected and 
experimentally measured number of events in the zth bin, while V denotes the covariance 
matrix. The statistical and systematic uncertainties for each bin are summed together in 
quadrature to form the diagonal entries af^ = a|(stat)-|-af (sys) in V, while the off-diagonal 
elements which take into account residual bin-to-bin correlations among the systematic 
errors are given by afj = cTi(sys)aj(sys) [|2^. We illustrate in fig. |] the dependence of 
for 22 degrees of freedom upon A~^ over the domain -0.35 TeV"^ < A'^ < 0.45 TeV"^ 
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using the MT set SL structure function. @ Various limits may readily be extracted from 
the clean parabola appearing in this plot. For instance, we find 

A-2 = 0.113 ±0.080 TeV"^ (19) 

by locating the Xmin + ^ points on the parabola. This translates into the asymmetrical la 
interval 

A = 2.98+Q^^ TeV (20) 

for the gluon compositeness scale. Alternatively, we may quote the more conservative lower 
bound 

A > 2.03 TeV at 95% CL. (21) 

Analogous results from the other leading and next-to-leading order distribution functions 
evaluated at the renormalization scales = and = are displayed in table 1. 

We see from the 95% lower limit entries in the last column of this table that the bound in 
( pTD represents a conservative estimate for A. It compares favorably with the CDF limit 
for the compositeness scale associated with quark substructure. 

It should soon be possible to substantially improve our limit on new gluon sector 
physics. The current 1992-93 Tevatron run is expected to collect a data sample five times 
larger than the one used in this analysis. Cross sections at higher transverse jet energies will 
be probed, and sensitivity to any nonrenormalizable operators in the effective Lagrangian 
will be enhanced. In the next few years, an integrated luminosity of 60 pb~^ is projected. 
We therefore look forward to updating our findings as new data comes forth from Batavia. 
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^ Negative values for A^^ do not imply imaginary values for A. Instead the phase of A^'^ is 
simply absorbed into the dimensionless C2 coefficient which multiplies the gluon operator O2 in 
our effective Lagrangian. 
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Distribution 
function 




normalization 
factor 


Amin 


A-2 / TeV-2 


A95/ TeV 


MT 


El/2 


1.35 ±0.01 


11.76 


0.113 ±0.080 


2.03 


set SL 


El 


1.61 ±0.01 


12.17 


0.081 ± 0.067 


2.29 


MT 


El/2 


1.17±0.01 


11.10 


0.104 ±0.085 


2.03 


set S 


El 


1.39 ±0.01 


11.34 


0.081 ± 0.071 


2.25 


MT 


El/2 


1.15 ±0.01 


11.15 


0.104 ±0.085 


2.03 


set Bl 


El 


1.37 ±0.010 


11.45 


0.081 ± 0.072 


2.24 


HMRS 


El/2 


1.08 ±0.01 


11.52 


0.081 ± 0.092 


2.08 


set B 


El 


1.28 ±0.01 


11.48 


0.077 ±0.077 


2.22 


CTEQ 


El/2 


1.25 ±0.01 


15.82 


-0.131 ±0.070 


2.02 


set L 


El 


1.48 ±0.01 


15.97 


-0.117± 0.067 


2.10 


CTEQ 


El/2 


1.21 ±0.01 


12.18 


0.014 ±0.084 


2.45 


set MS 


El 


1.43 ±0.01 


12.22 


0.018 ±0.071 


2.65 



Table 1 
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Figure Captions 



Fig. 1. Inclusive jet cross section plotted against transverse jet energy E^.. The data 
points are the experimental measurements reported by CDF. The solid curve 
represents the predictions of pure QCD with no composite interactions, while the 
dot-dashed curve illustrates the effect of gluon operator O2 with A = 2 TeV and 
C2(A) = —An. The theoretical results are based upon the leading order MT set 
SL distribution function evaluated at = E'^/2. 

Fig. 2. for 22 degrees of freedom plotted as a function of A~^. 
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